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The intention of this work is twofold, first to present a most simple system capable of simulating 
the intrinsic bosonic Josephson effect with photons, and second to study various outcomes deriving 
from inherent or external decoherence. A qubit induces an effective coupling between two externally 
pumped cavity modes. Without cavity losses and in the dispersive regime, intrinsic Josephson 
oscillations of photons between the two modes occurs. In this case, contrary to regular Markovian 
decoherence, the qubit purity shows a Gaussian decay and recurrence of its coherence. Due to 
intrinsic non-linearities, both the Josephson oscillations as well as the qubit properties display a rich 
collapse-revival structure, where, however, the complexity of the qubit evolution is in some sense 
stronger. The qubit as a meter of the photon dynamics is considered, and it is shown that qubit 
dephasing, originating for example from non-demolition measurements, results in an exponential 
destruction of the oscillations which manifests the coUectiveness of the Josephson effect. Non- 
selective qubit measurements, on the other hand, render a Zeno effect seen in a slowing down of the 
Josephson oscillations. Contrary to dephasing, cavity dissipation results in a Gaussian decay of the 
scaled Josephson oscillations. Finally, following Ponomarev et al. [Phys. Rev. Lett. 106, 010405 
(2011)] we analyze aspects of thermalization. In particular, despite similarities with the generic 
model studied by Ponomarev et al, our system does not seem to thermalize. 

PACS numbers: 03.65.Yz, 03.75.Lm, 05.30.Jp, 42.50.Pq 



I. INTRODUCTION 



Cavity quantum electrodynamics (QED) systems, and 
later circuit QED, operating in the strong coupling 
regime have proven to be excellent for studies of quan- 
tum phenomena such as entanglement/non- locality and 
non-classical states of the radiation field [H, [l]- In the 
recent past, however, new directions in these fields have 
been explored. A long term goal is to build many-body 
quantum models constructed of interacting photons. The 
idea is that these systems can simulate many-body mod- 
els appearing in fields such as condensed matter the- 
ory [3]. By now, several proposals of such cavity /circuit 
QED quantum simulators have been put forward, among 
others; Mott insulator-superfluid phase transitions [1], 
spin- models and fractional quantum Hall statesQ, the 
anomalous Hall effect @, relativistic effects 07 non- 
Abelian models Q, the Jahn- Teller effect @, and the 
photonic Josephson effect [HI- The desired non- 
linearity, and thereby photon-photon interaction, plays 
an important role in all of these and is established via 
interaction between light and matter. In the simplest ex- 
ample of matter-light interaction, the Jaynes-Cummings 
model (JC) [l^, the interaction is inherently non- 
linear giving rise to the effective scattering between pho- 
tons. This far, for the JC model the most striking outings 
of the non-linearity are the collapse-revival fl^ and the 
photon blockade effects [l^ • 

Experimentally, quantum simulators involving many- 



body systems are mainly to be found within ultra- 
cold atomic gases, either trapped dilute fermi or bose 
gases [l^ or more strongly correlated lattice or low di- 
mensional models [T3]- These systems provide an ex- 
treme controllability of system parameters as well as be- 
ing well isolated from their environments and hence pos- 
sessing long coherence times. Together with the experi- 
mental achievements within this area, topics such as de- 
coherence and thermalization in systems processing many 
degrees-of-freedom have gained extra attention recently. 
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FIG. 1: Schematic picture of the subsystem dynamics. The 
full A + B system evolves unitary under the map U (r). Trac- 
ing over subsystem B gives the state Pa{t) of subsystem A at 
the initial {t — 0) and final {t = r) times. The dynamical map 
L[pa,t] relates these two states of subsystem A. The evolu- 
tion of Pa{t) described by the dynamical map is in general 
not unitary, and even if pa{0) is pure the propagated state 
Pa('!") does not need to be pure. 



'Electronic address: 'j olarson@fysik.su. se| 



For completely isolated systems, it is understood that a 
true thermalization of the full system cannot occur for an 
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initial pure state. Nevertheless, expectation values (An), 
where An is any operator that operates in the subspace 
spanned by particle/subsystem n, often display charac- 
teristics as if the state of this particle/subsystem is 
a microcanonical one [l8|. Dividing the full system into 
two subsystems A and B, a sketch of the dynamics is 
given in Fig. [T] For more details see [Tol . [20j . The full 
system state p(0) evolves unitarily in time under some 
operator U{t). Tracing out the degrees of freedom of 
subsystem B at time t = or t = r gives the initial 
and final states for subsystem A, pa{0) and pa{t)- The 
state pa{0) is related to the state pa{t) via a dynam- 
ical map [ll, HOl- When the evolution of subsystem A 
is Markovian, the dynamical map L[pA, t] forms a semi- 
group [l^, llSl • The dynamics of subsystem A is in general 
not unitary. Moreover, when B contains a large number 
of degrees-of-freedom, any information in A is likely to be 
distributed over the many degrees-of-freedom of B, and 
PAiT) then becomes approximately diagonal in its com- 
putational basis. Regardless whether pAir) approaches 
a microcanonical state or some other diagonal state, sys- 
tem B has a decohering effect on A. If B is a closed 
finite system, we call this decohrence as inherent, while 
if the decoherence derives from coupling to an infinite 
reservoir we call it external. A general believe is that 
for closed quantum systems, chaos and non-integrability 
is deeply connected to thermalization, and thereby inte- 
grable many-body systems may not thermalize or even 
result in quasi diagonal states pa{t) In such sys- 

tems, recurrences of the coherence of system A often oc- 
curs, i.e. the information in system A, once lost to i?, is 
recovered in A at certain revival times. 

For the photonic many-body systems mentioned above, 
relatively little has been investigated in terms of inherent 
decoherence and thermalization. On the contrary, exter- 
nal decoherence and thermalization induced by a sur- 
rounding reservoir have been studied in great detail for 
cavity/circuit QED models [l^[22|. There are, however, 
some exceptions, both theoretical and experimental. In 
his seminal work (2^, J. Gea-Banacloche explained how 
inherent decoherence may cause peculiar effects in cav- 
ity QED systems. Considering the JC model with the 
field initially in a coherent state with a large amplitude, 
it was shown how the qubit disentangle from the field at 
half the revival time, regardless of the initial state of the 
qubit . In the qubit subspace, the evolution is clearly 
non- unitary since the overlap/angle between any two ini- 
tial qubit states is not preserved. A physical picture of 
this phenomenon is given in "245. In a couple of experi- 
ments performed by the ENS group in Paris, loss of qubit 
coherence due to interaction with a microwave cavity field 
was demonstrated [25]. A similar idea was later applied 
in circuit QED, and the qubit dephasing deriving from 
measurement backaction was studied [2^ . In the present 
work we analyze a plethora of decoherence phenomena, 
both inherent and external, in a bimodal cavity /circuit 
QED system. Decoherence originating from; non-unitary 
subsystem evolution, irreversible coupling to surrounding 



environments, as well as irreversible coupling to a mea- 
surement device are discussed. 

In particular, a two-level system (qubit) induces an ef- 
fective coupling between two quantized photon modes. 
In the dispersive regime, as is studied here, the qubit 
mediated tunneling between the two boson modes shows 
much resemblance with a bosonic Josephson junction, 
appearing for example in Bose-Einstein condensates in 
double- well potentials [13, HI] or via coupling of atomic 
Zeeman levels in a spinor condensate [l^. The photonic 
counterpart of the Josephson effect has currently been 
considered in either coupled cavity systems (where the 
non-linearity is of a JC photon blockade type) [l^l or 
in optomechanical cavity systems (where non-linearity 
instead stems from a mechanical photon blockade ef- 
fect) Unlike Refs. [13, the Josephson effect 
studied in this work is intrinsic and not as in earlier pub- 
lications where the effect is extrinsic and the Josephson 
current appears in real space. Moreover, the qubit in 
our scheme brings out both the boson tunneling as well 
as the non-linearity. In this respect, our model shares 
more similarities with the recent works of Ref. 30] . The 
models of [sQ are, however, both linear. Studying the 
linear entropy or the qubit Bloch vector, we find that 
inherent decoherence shows a Gaussian decay signaling 
non-Markovianity. The non-linearity causes a collapse- 
revival pattern in the Josephson oscillations, similar to 
the collapse-revivals found in the JC model [12] ■ In con- 
trast, the linear entropy of the qubit shows a much richer 
collapse-revival structure containing a series of well es- 
tablished fractional revivals. Measuring the qubit state, 
and especially the Bloch vector amplitude, provides di- 
rect information about the Josephson dynamics. How- 
ever, the qubit dephasing resulting from non-demolition 
qubit measurement causes an exponential decay in the 
oscillations. For non-selective qubit measurements, the 
Josephson current can be considerably slowed down via 
the quantum Zeno effect. Somewhat surprising, for a 
dissipating cavity, the scaled Josephson current does not 
decay exponentially, but rather Gaussian and it is found 
that for realistic system parameters around 10 full os- 
cillations survive. The bimodal JC model is integrable 
and thermalization is not to be expected. Including the 
counter rotating terms, i.e. going beyond the rotating 
wave approximation, we consider a weaker sort of ther- 
malization related to the analysis of [3l|. Despite the 
close similarity with the generic model of Ref. |31], wc do 
not encounter thermalization. We further evidence the 
absence of thermalization by consider long-time state av- 
erages as was recently proposed in [s^]- Neither in this 
sense does our system thermalize. 

The paper continues as follows. In the next section, the 
model system is introduced and the photonic Josephson 
effect demonstrated. In the linear regime we present an 
analytical approximation giving a deeper understanding 
of the effect. Section Hill proceeds by considering the evo- 
lution of the qubit subsystem. Since the full system is so 
far assumed closed, the field and the qubit subsystems ex- 
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hibits the same type of reduced evohition. Gaussian de- 
cay together with periodic recurrences of the coherence 
is found. Effects deriving from non-demohtion or non- 
selective qubit measurements are addressed in Sec. IIV[ 
while the following Sec. |V] considers the more complete 
picture where the cavity has been coupled to a zero tem- 
perature photon bath. In Sec. I VII the prospects of a 
weak kind of thermalization is analyzed. Finally, con- 
cluding remarks are presented in Sec. IVIIl 



II. INTRINSIC PHOTONIC JOSEPHSON 
EFFECT IN A BIMODAL RESONATOR 

We consider a two- level system, either an atom/ion in- 
teracting with an optical Fabry-Perot cavity or a super- 
conducting qubit embedded in a transmission line res- 
onator. We will call the two-level system a qubit and 
describe it using the Pauli matrices [(Tz,i7''=] = ±2g^ 
and a^lib) — ±|±), where |±) are the two internal qubit 
states. The energy difference between the two qubit 
states is Tifi. The resonator supports two quantized 
modes, quasi-resonant with the qubit transition. We will 
take the two modes as degenerate with frequency Wc, but 
point out that the degeneracy is not crucial and could in 
principle be lifted. Introducing coherent driving of the 
two modes with amplitude 77, the Hamiltonian becomes 



H = uj. 
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irj 
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(1) 

where d and b (a^^ and 6+) are the annihilation (creation) 
operators of the two modes, ga,b are the qubit-field cou- 
pling strengths that have been taken real, and uip is the 
frequency of the driving field. The third line represents 
the counter-rotating terms. Turning to a rotating frame 
with respect to the pumping, and after imposing the 
rotating- wave approximation (RWA), we have the time- 
independent Hamiltonian 



H. 



RWA 



= S (^d^d + b^bj + y(7^ 

+5a (a+CT" + a+d) + gt (h^'a^ + a+bj 



(2) 
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Here, S = ^ and A ~ fl — ujp are the cavity- 
pump and atom-pump detunings respectively. For res- 
onant driving {5 = 0) of a cavity without the qubit. 



energy will be constantly pumped into the two modes. 
For (5 ^ 0, on the other hand, the system will reach a 
coherent steady state with {fia.b) — v'^/S'^. 

With g = y/g^^+gl and cos0 ~ ga/g, we define a new 
set of boson operators 
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For ga — gt, which will be the situation considered 
throughout this work, the above rotation is a simple 
Hadamard rotation of the two boson modes. In the loss- 
less case, assuming zero pumping (77 = 0), the RWA 
Hamiltonian is 



Ho = S 



(iti- 



B^B 



+gV2(A'<a- +a+A 



(4) 



Thus, in the transformed picture the two boson modes 
are decoupled and the qubit interacts with the A-mode 
according to a JC type. In the large detuning regime. 



Ho = S (A'<A+B^B] + 2UA'iAa 



A ^ g\J{A''A), we impose an adiabatic approxima- 
tion [33I I to find, up to some trivial constants, 

A 

(5) 

where U = g^/A. Even though the Hamiltonian is diago- 
nal in the Fock basis of the A- and B-modes, the original 
a- and 6-modes still couple in a nontrivial manner since 



A^A ^ ^ + fib + d^b + 6+a) , 



(6) 



where the photon numbers fia = a+a. Neglecting the 
non-linear terms, the linearized Hamiltonian becomes 



^0 ~ {fia +nb) + U (^d^b + fetfl) 



(7) 



with the Stark shifted frequencies 5^ — S + azU. The 
total number of photons, N = fia + fib, is conserved and 
the dynamics depends solely on U. 

With the Hamiltonian Hq of Eq. ([7]), it is clear that 
any initial field state \^pa,0), where {ipa, V'b) represents the 
composite state with the a-mode in \ipa} and the 6-mode 
in \ipb), transforms into (up to an overall phase) 10,-00) 
after a time tj — tt/2U, regardless of the atomic ini- 
tial state. This is the analog of intrinsic bosonic Joseph- 
son oscillations for photons. In the earlier works on the 
photonic Josephson effect [13, [lil photons tunneled be- 
tween spatially separated regions, while here they tunnel 
between two orthogonal cavity modes. To make things 
more transparent, lets assume the a-mode to be initially 
in a coherent state with amplitude a, d\a) = a\a), the 
6-mode in vacuum, and the qubit in the superposition 
state \ip) = (1+) + |-))/\/2, i.e. 
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FIG. 2: Time evolution of the scaled photon inversion Z(t) 
(a) and (c), and linear entropy (b) for the Hamiltonian ((2|. 
As an initial state we have a coherent state with amplitude 
Of = 4 in the a-mode, vacuum in the 6-mode, and the qubit 
in (1+) + \ Due to the non-linearity of the JC model, 

a clear collapse-revival pattern emerges. The linear entropy 
displays a more pronounced collapse-revival pattern than the 
inversion. The dimensionless parameters are jr = 1, A = 50, 
5=1, and r] = 0.1. 



propagated state return back in phase a revival of the 
Josephson oscillations occurs. For revivals to be possi- 
ble, the spectrum should be discrete and quasi-linear, 
i.e. the non-linearity cannot become to strong. In the 
language of cavity QED, revivals is a direct pro of of the 
quantized nature of the electromagnetic field [33] , while a 
collapse could, in principle, derive from any fluctuations. 
We conclude from the fact that our system is governed by 
JC dynamics, a self-trapping effect should not be found, 
while collapse-revivals are expected. 

The result for the scaled photon inversion Z{t), ob- 
tained from numerically solving the driven two-mode JC 
model (12), is presented in Fig. [21 We truncate the Hilbert 
space such that Tr[p(t)] w 1, with p{t) the full density 
operator. In (a) we give the inversion for long times, 
while in (c) we display the short time behavior. The fast 
oscillating fluctuations in Z[t) seen in (c) derives from 
the non-zero driving (77 ^ 0). The long time dynamics 
exhibits a typical collapse-revival structure as found in 
the JC model [l^, [3, IS]. However, in the present bi- 
modal model the collapse-revival structure is much more 
complex compared to the regular JC model. In order to 
demonstrate this we introduce the linear entropy of the 
field or the qubit [11] 
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In the interaction picture, the time evolved state reads 
1 
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and for the scaled photon inversion defined as 

(ha) + (hb) 



(9) 



(10) 



we have Z{t) — cos{2Ut). This result is an ideal case 
where losses and non-linear contributions have been ig- 
nored. In the next sections we will discuss the influences 
of dissipation and decoherence, while we now consider 
the consequences of the non-linearity. 

Typically, non-linearity plays an important role in 
Josephson dynamics. In the bosonic Josephson junc- 
tion, it gives rise to effects like self-trapping [l^, [13 and 
collapse-revivals [SJ, |35| . Self-trapping appears since the 
non-linearity can be seen as inducing effective shifts in 
the two mode frequencies and when these are out of res- 
onance the population swapping is hindered. A collapse 
in the oscillations derives from the fact that the non- 
linearity disturbs the equividistant energy spectrum of 
the linear model ([7]) and this in turn leads to a destruc- 
tive interference effect in the dynamics [36j . If at some 
later time t^, the evolving eigenphases contained in the 



where pf^q = TTqj[p] is the reduced density operator of 
the field "/" or the qubit "g" and p is the full system's 
density operator. Linear entropy is a measure of pu- 
rity of the state p/^^, with Pf ,q = representing the pure 
state, and for the special case of a qubit maximum mixed- 
ness gives Pq — 1/2. Whenever the evolution is uni tary 
and the initial state is pure one has Pf{t) — Pq(t) [40]. 
Therefore, we do not need to distinguish between the 
two subsystem's linear entropies in this case. For the 
above example, we plot in Fig. [5] (b) the corresponding 
linear entropy. This time the collapse-revival structure 
is much more complex, indicating that rich dynamics is 
occurring even within the collapse region indicated by 
(a) . The revival-peaks are marked according to standard 
conventions [s^. In a collapse region, the phase-space 
distributions are smeared out over the whole accessible 
part of the phase-space (ilj. At the 1/2 peak the distri- 
bution reforms, however 180 degrees out of phase from 
its initial position, while at the 1/1 revival it reforms 
at its original location [4l|. The remaining fractional 
revivals represent partial restoration of the phase space 
distributions, which then form several localized blobs in 
phase space and not single ones as is the case for the 1/2 
and the 1/1 revivals [41|. Increasing the non-linearity, 
i.e. increasing the fraction g/A, results in shorter revival 
periods but also in less pronounced revivals. 
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III. FIELD INDUCED DECOHERENCE OF THE 
QUBIT 

Decoherence of a quantum system is seen as loss of in- 
formation, which is only possible by coupling the system 
to some other system. If the additional system contains 
a large number of degrees-of-freedom, the information 
can be shared in many ways and it takes a long time 
to restore the information in the initial system. If the 
number of degrees-of-freedom goes to infinity and the 
various coupling strengths varies, we expect the infor- 
mation to be forever lost. When we neglect any back- 
action of the reservoir onto the system we say that the 
evolution of our system is Markovian [l^. Markovian 
dynamics is described by irreversible evolution (i.e. non- 
unitary), typically leading to exponential loss rates. As 
the reservoir size is finite or for very structured reser- 
voirs, non-Markovian effects may become important giv- 
ing rise to recurrences of information, and moreover to 
non-exponential decay [l^. In particular, for short times 
a Markovian decay is linear in time while non-Markovian 
decay is typically quadratic. In this section we will dis- 
cuss the inherent decoherence effects the cavity fields in- 
duces on the qubit. 

In some recent publications, Sokolovski and co-workers 
considered a quibit coupled to a bosonic Josephson junc- 
tion [s^l- A Bose-Einstein condensate trapped in a 
double- well potential is analyzed. The qubit couples to 
the barrier separating the two wells, and in this respect 
its internal state controls the tunneling rate between the 
condensates. As such, they introduced the name "quan- 
tum gatekeeper model" for this setup. At first, the con- 
densate was seen as a reservoir for the qubit, and the in- 
herent decoherence of the qubit was studied. The qubit 
as a measuring device for the bosonic Josephson oscil- 
lations was then analyzed. In many aspects, the cur- 
rent system shows great resemblance with the gatekeeper 
model. There are some crucial differences though; the 
bosonic modes in our model induces coupling between the 
qubit levels, i.e. our model is intrinsically non-linear, and 
the photon tunneling is present for both internal qubit 
states. In principle, as demonstrated in the previous sec- 
tion, transitions within the qubit can be neglected by 
considering the limit of large detuning. Moreover, by in- 
cluding a second qubit residing all the time in its lower 
state |— ) but still couple identically to the boson modes 
as the first qubit, it is possible to construct an identical 
gatekeeper model as in [s^. We, however, will keep only 
a single qubit and allow for qubit transitions and thereby 
for non-linear effects. 

Throughout this work, the two qubit states |±) are 
assumed meta-stable such that dissipation and decoher- 
ence can be neglected. In this section we further take 
any cavity losses to be zero within the time-scale of the 
experiment. Cavity losses and its effects will be investi- 
gated in Sec. IVII Our initial state will be taken to have 
the form ©. 

The time evolution in the dispersive linear regime is 
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FIG. 3: (Color online) Comparison of the linear entropy as 
obtained from the analytical approximate result ((14} or from 
numerical diagonalizition of the Hamiltonian ((2} (analytic ap- 
proximation (green lines) and numerical (black curves)). In 
(a) the parameters are as in Fig. [2l i.e. the non-linearities 
are weak, while in (b) A = 20 increasing the effects of non- 
linearities. 



given by Eq. ([9]), we find the reduced density operator 
for the qubit 



P++ P-+ 
p+- P-- 



where 



= 1/2 and 



P-+ 



lg-|a|^sin^(2;7t)gi2|a|^ sin(2(7t) 



such that the linear entropy becomes 
Pit) = 1 (l - e-^l" 



^{2Ut) 



(12) 



(13) 



(14) 



The short time decay of the coherence term p ^ is 

quadratic in t as expected for non-Markovian decoher- 
ence, and similarly the linear entropy grows quadratically 
in t. In this linear approximate model, at i = 7r/2[/, 
the coherence is fully recovered. For unitary evolution 
with initial pure states, the linear entropy is a good 
measure of entanglement, and we therefore see that the 
two fields quickly become entangled with the qubit. For 
large field amplitudes |a|, the two subsystems stays al- 
most completely entangled until they fully disentangle at 
tj — 7r/2[/. Thus, the qubit rapidly loses information 
(coherence) to the two boson modes, but since the sys- 
tem is finite the information is restored periodically in 
time. Once the information is completely transfered to 
the two modes, the qubit is maximally entangled with the 
cavity fields. Note that the disentanglement occurs when 
the photonic Josephson junction has performed half an 
oscillation, i.e. when the initial coherent state in mode a 
has been transfered to mode 6. 

The above analytical model is an idealized case where 
we completely disregarded non-linearities. In Fig. |31 we 
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therefor compare the analytical results with results ob- 
tained from a full numerical diagonalization of the Hamil- 
tonian. As for Fig.[3l the non-zero pumping induces small 
rapid oscillations. In the upper plot (a), the parameters 
are the same as in Fig. [2] and we find that the analytic 
approximations captures well the overall structures. The 
oscillations deriving from the pumping is naturally not 
seen in the analytic results. In the lower plot (b) we 
decrease the detuning A in order to increase the non- 
linearity, and as expected the two results begin to dis- 
agree. Figure [3] gives dynamics at the time-scale of the 
Josephson oscillations, while in the previous section we 
studied evolution on the scale of collapse-revivals. Nat- 
urally, the linear analytical model of Eq. (|14|) is not ca- 
pable of catching the collapse-revival patterns of Fig. [2] 
(b). 




FIG. 4: (Color online) The integration of the Bloch vector 
length obtained from the numerical results (black curve) and 
from the analytical expression (I16p (green curve). The param- 
eters are taken the same as in Fig. [3] (a). When integrated, 
the discrepancy between the analytical and numerical results 
for the Bloch vector amplitude adds up and the two curves 
disagree, even though the qualitative structure is the same. 



vector. In our analytic model we have 

((7,) = e-l"l'^'"'(2C/t) cos[2|a|2sin(2?7i)], 

(15) 

{cjy) = e-l"l'^'"'(2C/t) sin[2|a|2sin(2C/t)] 

and for the length of the Bloch vector 

|i?(i)| =e-l°l'-"'(2t^*). (16) 

In the vicinity of the disentanglement times tj = mr/2U 
{n = 0, 1, 2, ...), the Bloch vector length is roughly 
unity. In between the disentanglement, the vector length 
rapidly shrinks to approximately zero. The area cov- 
ered by the Bloch vector during one disentanglement pe- 
riod can be approximated by / = J die"'"' ^^^^^ = 

^7r/4[/|a|2. The function e{t) = /J dt'\R{t')\ should 
thereby have a stair-like shape increasing in steps ~ /. In 
Fig.|3]we present the numerical results of Q{t) compared 
to the integration of \R{t)\ of Eq. (|16p . The parame- 
ters and initial conditions are the same as for Fig. [3] (a) . 
From the figure it is clear that even in the collapse re- 
gions, the two subsystems (fields and qubit) are not fully 
entangled, which would be reflected in a Bloch vector 
of zero length. The analytic results (green curve) still 
captures the general structure of the numerical results 
(black curve). As a conclusion, by measuring |i?(t)|, and 
thereby Q{t), the number of Josephson oscillations per- 
formed by the cavity fields is obtained. Measuring the 
qubit dipole movements (ax) and (ay) is easily obtained 
by detecting the qubit population inversion after a 7r/2 
pulse has been applied to the qubit [1]. Any measure- 
ment on the qubit will however induce a perturbance, 
whether it is in the qubit population or in its coherence. 
While such effects were disregarded from the analysis of 
Ref [S^l, we address them in the next section. 

IV. QUBIT MEASUREMENT INDUCED 
DECOHERENC OF THE JOSEPHSON EFFECT 



In the "gatekeeper model" of [s^l, the qubit was sug- 
gested as a quantum "meter" measuring the Josephson 
oscillations of the Bose-Einstein condensate. In that 
model, if the qubit occupies its upper state |-|-) Joseph- 
son oscillations took place, while if the qubit resided in 
its lower state |— ) oscillations were blocked. Thereby, 
by measuring the effective amount of time the qubit has 
spent in the state |-|-) one knows for how long the oscil- 
lations ("the gate") has been turned on ("open"). The 
situation is different in our case since Josephson oscil- 
lations are independent of whether the qubit is in |+) 
or |— ). This, however, does not imply that the qubit 
does not carry any information about the oscillations. 
In a Bloch picture, since (ct^) « the Bloch vector 
R(t) = (((Tx), (o'jf), (o'z)) revolve mainly along the equa- 
tor of the Bloch sphere. Thereby, in our system it should 
be more practical to measure (ax) or (ay) instead of ((T^). 
Another possibility is to measure the length of the Bloch 



In its original idea 142] , the Josephson effect is a macro- 
scopic quantum phenomenon describing a current across 
the bulk of two attached superconductors. A common 
quantum coherence among the constituent electrons is 
the crucial aspect for a non-zero current. At the mean- 
field level, the relative phase between the two supercon- 
ductors are directly coupled to the inversion Z{t) 
The question in our system then raises how a qubit mea- 
surement will affect the Josephson oscillations. For ex- 
ample, we note that since the Stark-shifts of the qubit 
states |±) depends on the photon number, the back- 
action noise from measuring the relative phase between 
the qubit states should manifest itself in the phase co- 
herence of the cavity fields. Apart from dephasing, more 
drastic outcomes may appear when the measurement is 
projective. Most well known is the Zeno slowing down 
caused by repeatedly reprojecting the system state. In 
this section, we will analyze both the case of pure de- 
phasing and Zeno slowing down. 
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A. Measurement induced dephasing 

A qubit measurement relying on for example ioniza- 
tion is, of course, destructive. Once measured, the qubit 
is "destroyed". The idea of a quantum non-demolition 
(QND) measurement is that repeated measurements can 
be performed on one and the same system, all yielding 
the same outcome. Let us take a physical quantity char- 
acterized by some operator A, then if A commutes with 
the system Hamiltonian and the measurement of A does 
not induce any back-action noise on A, a, second measure- 
ment of A will reveal the same value as the first measure- 
ment 43]. Even though the measurement does not have 
any back-action effect on A, according to quantum me- 
chanics it will induce noise in the conjugate variable of 
A. In cavity QED, in order to measure the photon num- 
ber of the quantized cavity mode, QND measurements 
has been thoroughly analyzed both theoretically [31 and 
experimentally |45| . Similar QND measurements have 
been carried out also in circuit QED [1^1 ■ Our goal is 
not, however, to measure the number of photons con- 
tained in the resonator, but to measure the state of the 
qubit. For superconducting qubits, such QND measure- 
ments have recently been demonstrated 47| . The picture 
we have in mind is that the the qubit level-population 
is detected non-destructively, i.e. measuring ct^. Since 
neither ax nor cTy commutes with ct^, the measurement 
induces a qubit dephasing [i^. In the bosonic Joseph- 
son junction, decoherence has been shown to suppress 
the oscillations [H, . For Bose-Einstein condensates 
in double- well potentials, the loss of coherence (phase 
diffusion) is mainly due to atom-atom interactions [50| . 
The dephasing we consider in this section has a different 
nature as it affects the anchilla system and not the boson 
modes directly. 
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FIG. 5: (Color online) The scaled photonic inversion Z{t) 
of Eq. (|10p and the amplitude of the Bloch vector |-R(t)|, as 
obtained from numerical integration of the master equation 
(|17|l (black curves). The green curve shows the comparison 
with the analytical approximation ((TSJ for the Bloch vector 
amplitude. Both the Josephson oscillations and the Bloch 
vector decay exponential due to the qubit measurements. The 
decoherence rate 7 — 0.005, and the remaining parameters are 
as in Fig. [2l 



We assume this QND measurement to be an irre- 
versible process that can be modeled by a master equa- 
tion [1^, i.e. the measurement is continuous and far from 
a traditional von Neumann measurement. As no dissipa- 
tion occurs, we write down the equation-of-motion for 
the full density operator as 



dpjt) 
dt 



~i[p{t),HRWA] 



(17) 



xl ' 



where 7 is the effective decoherence rate induced by the 
measurements. Letting 7 = we find identical evolution 
as in Figs. [2]and|3l Once 7 7^ 0, the off-diagonal terms 

p |_ and of the atomic density operator will typically 

show an exponential decay [l9| . 

We solve the master equation ([T7)) by simple New- 
ton integration, valid to order 0{dfi) in the time-step 
dt. However, in its simplest form, the first part describ- 
ing regular Hamiltonian dynamics becomes non-unitary 
and the method is greatly unstable. Therefore, the first 
unitary evolution is taken into account in all orders of 
dt^ while the irreversible evolution is performed to order 
0{dt^). As long as dt is small this should be justified, 
and to be fully sure we have verified convergence of our 
results by varying dt. 

In Fig. m the photonic inversion (a) together with the 
Bloch vector amplitude (b) are displayed. It is seen that 
the Josephson oscillations decay exponential, indicating 
that the noise induced by the qubit measurement propa- 
gates to the fields and thereby demolish the macroscopic 
coherence needed for the oscillations. As pointed out in 
the introduction, this exponential decay signals an irre- 
versible loss of coherence. The Bloch vector amplitude 
shows an interplay between the non-Markovian inherent 
decoherence resulting from coupling to the two cavity 
modes and Markovian external decoherence due to the 
qubit measurements. The exponential decay rate is set 
by 7, and in the most simple approximation we get 



\R{t)\ 



sin^(2;7t)-7t 



(18) 



From Fig. [5] (b) we note that the analytical approxima- 
tion ([T5|) (green curve) slightly overestimates the expo- 
nential decay. A large 7 implies that information about 
the qubit, and thereby the boson modes, is rapidly leaked 
from the system to the "observer/meter" . If the mea- 
surement is perfect, i.e. a complete knowledge of the 
system is obtained, all the coherence shared between the 
two subsystems gets lost, and the qubit disentangles from 
the two modes. There is hence a trade-off between the 
information subtracted and the survival of the Joseph- 
son oscillations. We note that using the qubit as a meter 
was recently utilized in a circuit scheme measuring the 
quantum state of a mechanical (rather than an optical) 
resonator [51 1. 



8 



B. Quantum Zeno effect 

Repeatedly measuring a system affects the effective 
evolution of the system. The extreme situation of freez- 
ing the evolution by constantly projecting the system 
state back onto its initial state is the idea behind the 
quantum Zeno effect (s^ . In reality, no measurement is 
expected to be instantaneous resulting in a "dead-time" 
after each measurement. In this respect an ideal Zeno 
effect is unlikely f52!|. This, however, does not prevent 
a slowing down effect, which was recently demonstrated 
within cavity QED |5^. We briefly discussed the evolu- 
tion of a cavity mode under coherent external pumping. 
For an empty mode, the field growth is initially quadratic 
in time. Under a Zeno slowing down, this becomes linear 
in time [s^. A quadratic growth hinges on coherently 
adding up the injected field, while repeated QND mea- 
surements destroys the field phase. The injected field can 
then be seen as making a random walk in phase space and 
thereby the linear growth. We will not consider the field 
build up due to external pumping, but instead consider 
the growth of the initially empty cavity mode deriving 
from photon tunneling between the two modes. We have 
seen that such growth is as well quadratic for small times 
t. In order to separate it from external pumping we pick 
ry = 0. 

The approach of the previous subsection is not capa- 
ble of describing a Zeno slowing down since we pictured 
a constant loss of information, and not repeated phase 
disturbances. We will therefore instead assume repeated 
non-selective qubit measurements [13] ■ At times t„ = nr 
(n = 1, 2, ...) the qubit population is measured but the 
result is not recorded. The correlations between the fields 
and the qubit is taken to be fully demolished by the mea- 
surement, hence 

p{t) ^ Pf{t)p,{t). (19) 

In between measurements, the system evolves unitarily 
as Ur = exp (^—iHfiwATj . Since the measurements are 
taken to be instantaneous, the full evolution operator 
becomes 

r 1 k 

U{t, t) = MUr , (20) 

where M is the field-qubit disentanglement operator 
characterizing the measurement, and k is the largest in- 
teger such that < t — /cT < 1 . If r is small, the coherent 
build up of the field is frequently interrupted and we ex- 
pect a slow and linear field increase. 

The predicted linear rise of the field amplitude is nu- 
merically demonstrated in Fig. [5] for various values of r. 
For no measurements, r = oo, the field grows quadrat- 
ically from time t — 0. When r is decreased we find a 
regime of linear growth which passes over into a quadratic 
regime. The smaller r, the smaller is the incline and the 
more inhibited is the quadratic evolution. The figure is 
a typical example of a Zeno slowing down [s^ [sH] • 




FIG. 6: (Color online) Evolution of the field intensity {hi,} un- 
der repeated non-selective qubit measurements. The different 
curves correspond to different measurement frequencies; red 
curve T — oo, green curve r — 1/100, blue curve r — 1/500, 
pink curve r = 1/1000, and black curve r = 1/1500. The 
dimensionless parameters are, g = 1, A = 50, 5 — 1, r] — 0, 
and Q = 2. 



For a composite system A + B, it has been demon- 
strated that repeated projective measurements on sub- 
system A render a purification of subsystem B f56| . 
Letting the p(0) = pa(0)pb(0), with ^^(0) = 
\4>a{0)){'4'a{Q)\, be the initial state, the projective mea- 
surements divides the unitary evolution at tn = nr 
(n = 1, 2, ...) according to 



p(t)->PB(t)|V'A(0))(V^A(0)|. 



(21) 



In the case of non-selective measurements a disentangle- 
ment between the two subsystems appears, while for pro- 
jective measurements in addition to disentanglement the 
A system is as well reset to its initial state. Regardless of 
the initial state pb{0), the effective evolution will trans- 
late into a pure state. This is the idea behind measure- 
ment controlled evolution [s^ Is^j ■ In the present system, 
we have verified that this mechanism indeed works. Var- 
ious initial states like coherent and thermal field states 
have been considered, and in every case the purity of the 
field subsystem approaches unity after sufficiently many 
projective qubit measurements. The specific purified field 
state differs between various initial field states. In cer- 
tain situations the fields approach Fock states, but this 
seems not to be a generic feature. 



V. RESERVOIR INDUCED DECOHERENCE 

So far we have considered decoherence as either inher- 
ent or imposed via qubit measurements. We now proceed 
by analyzing the effects deriving from coupling the pho- 
ton modes to a zero temperature reservoir. The two qubit 
states |±) are still taken to be metastable and losses are 
therefore solely in the boson fields. As for the case of 
qubit measurements, we model the irreversible coupling 
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by a master equation [22 

dpit) 



dt 



+K (2ap(t)at - _ p{t)na) (22) 



Here, k is the common decay rate of the two modes. The 
above equation is commonly used for a zero temperature 
photon reservoir, the generahzation to non-zero temper- 
ature typically implies two extra Lindblad terms that are 
as well proportional to k. In the case of strong qubit-field 
coupling, it has been shown that such a master equation 
does not in general decay into the correct canonical distri- 
bution [ssj . However, at zero temperature as considered 
here, the long time behavior of the above master equation 
is physically relevant [58.]. In contrast to Eq. p7|) which 
describes decoherence, the above equation ([^ results in 
dissipation of cavity photons to the reservoirs. As will be 
seen, pure dephasing and dissipation have different influ- 
ences on the Josephson oscillations. Direct dephasing of 
the boson modes, without dissipation, is obtained by re- 
placing a and by ha (and similarly for the b mode) in 
the two Linblad terms. 
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FIG. 7: Evolution of the scaled photon inversion Z{t) (a) 
and the Bloch vector amplitude |-R(t)| (b) in the presence of 
photon dissipation. The Bloch vector amplitude displays a 
combination of Gaussian and exponential decay. The scaled 
Josephson oscillations, however, survives much longer. The 
parameters are as in Fig. [2l but with a = 3 and k = 0.0012(; 
corresponding to the experimental decay rate of Ref. [goI ]. 



Exactly like in the previous sections, we can find an 
analytical solution of Eq. (|22l) in the linear dispersive 
regime [59i]. In the long time limit and letting ^ = 0, the 
atomic density operator becomes 



Pat[t) = - 



1 P+- 
P-+ 1 



where 



P+- 



(23) 



(24) 



and p is a complex time-independent factor of order 
unity [1^ . The form of the field density operator reads 

Pf = l [\a+,b+){a+,b+\ + |a_,&)(a-,6-|] , (25) 

where the explicit expressions for the coherent ampli- 
tudes a± (t) and b + (t) in the most general case are rather 
complex, see Ref. [59|. They become, however, simple in 
the limit of no pumping, ry = 0, 



W = Y [l + cos(2[/t)]. 



= — e-2«* [l-cos(2C/i)]. 



(26) 



It follows that the inversion Z{t) = cos{2Ut) as without 
losses. In general though, with rj ^ 0, the expression for 
the inversion becomes very complicated. Nonetheless, for 
this case we note that generally {fia) 7^ (hb) in the limit 
t 00. This derives from the fact that the two mode 
were initially differently populated. As a consequence 
of this, the photonic inversion does not approach zero; 
limt_j.oo Z{t) ^ 0. Moreover, as long as (5 7^ 0, the initial 
decay of Z{t) is not exponential, which is in contrast 
to pure dephasing. We will also see that non-linearities 
beyond the above analytical model play an important 
role. 

As before, our numerical investigation goes beyond the 
linear regime of the analytical consideration. We solve 
the master equation (|22p with the same procedure as in 
the previous section. For the photon decay rate we take 
an experimentally relevant value of k = 0.0012(7 The 
numerical results are shown in Fig. [71 The rest of the pa- 
rameters are the same as in the example of Fig. [51 except 
that the amplitude of the initial coherent state a = 3 
(instead of a = 4 as in Fig. [2]). As in the case of pure 
dephasing, the Bloch vector amplitude displays an inter- 
play between Gaussian decay and recurrences and expo- 
nential decay. Contrary to the example of dephasing, the 
amplitude does not approach zero in the long time limit. 
This does not, however, imply that the qubit is entangled 
with the boson modes. Since the evolution is not unitary, 
the linear entropy (nor the Bloch vector amplitude) is a 
good measure of entanglement. In fact, the off-diagonal 
elements of both pat and pf vanish in the long time limit. 

It is interesting to note that the behavior of the sys- 
tem with photon dissipation is very different from the 
case with dephasing. Dephasing induced an exponential 
decay of the photon numbers (ha) and {fib) such that 
(fia) = (hb) for large times. The same kind of exponential 
decay was seen in the photon inversion Z{t) (see Fig. [5l 
(a)). Dissipation causes as well an exponential decay of 
(ha) and (hb) towards their steady state values, but Z{t), 
on the other hand, does not decay exponentially. At first, 
the scaled Josephson current is not greatly affected by 
the dissipation. After ^ 6 full Josephson oscillations, the 
scaled current begin to decay, however with something 
more reminiscent of Gaussian decay. As the amplitudes 
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of the two modes decrease, the rapid oscihations deriving 
from the non-zero pumping becomes more pronounced, 
and after around 20 full oscillations the pumping induced 
fluctuations washes out the Josephson oscillations and a 
new structure in Z{t) emerges. This additional inference 
pattern is not seen in the regime where the Josephson os- 
cillations are still well resolved. It is important to point 
out, even though the Josephson oscillations may survive 
relatively long, it should be remembered that the photon 
amplitudes {fia) and (fif,) can become very small putting 
limitations on photon detection. Today, on the other 
hand, very small field amplitudes ((n) ^1) can be mea- 
sured ^li]. 

In order to confirm the differences between dephasing 
and dissipation of the boson modes (Figs [5] and [7]) , we 
have numerically solved the dephasing master equation, 
i.e. replacing a and by fia and replacing b and &^ by 
fii, in the Linblad terms. This time, the inversion Z{t) 
decays exponentially in a similar way as in Fig. [5] (a). 

Figure [7] has been calculated using physically relevant 
parameters and we thereby conclude that the intrinsic 
photonic Josephson effect should indeed be realizable. 
However, only around 10 Josephson oscillations survive. 
On the other hand, also for the intrinsic and extrinsic 
bosonic Josephson effect studied in Bose-Einstein con- 
densates is typically only detected during a few oscilla- 
tions [1^ [6l| . Our numerical analysis does not allow us 
to go to larger photon numbers, but increasing the initial 
amplitude a may prolong the oscillation time. It should 
be noted though, that increasing the number of photons 
while keeping the detuning A fixed increases the non- 
linearity, which might lead to undesired effects. Stronger 
pumping will help in terms of preserving a measurable 
number of photons in the two modes, but it will as well 
imply a bigger impact of rapid oscillations on top of the 
Josephson oscillations. Optimizing the system parame- 
ters might be tricky, but it is nonetheless clear that the 
example of Fig. [7] is probably not at all the most optimal 
one in terms of a sustainable Josephson effect. 

VI. QUANTUM THERMALIZATION? 

Thermalization in quantum systems is normally at- 
tributed to dissipation and external decoherence [l^, HI] . 
Assuming a quantum system S to be irreversibly coupled 
to a thermal reservoir with some characteristic temper- 
ature T, and denoting the state of the system as p^, by 
waiting long enough the state of the system normally be- 
comes diagonal and on the form 

PsO,e-^^% (27) 

where Hs is the system Hamiltonian and /3 — fc^T is 
the Boltzmann's constant time temperature T. The re- 
duction of the off-diagonal coherence terms of ps is an 
outcome of the dissipation and decoherence induced by 
the reservoir. On the other hand, assuming ps to be ini- 
tially pure and the system to be closed, it follows that 



Ps will remain pure, i.e. Tr[/9^(t)] = 1. However, this 
does not forbid the state pA of subsystem A of the to- 
tal system S to become diagonal. We saw examples of 
temporary subsystem diagonalization in the previous sec- 
tions. This, of course, is an effect of the inherent decoher- 
ence enforced upon subsystem A. Whether a recurrence 
of the off-diagonal terms appears is a highly legitimate 
question. More precisely, one wishes to understand if, 
why, and how pA attains an asymptotic diagonal distri- 
bution that can be modeled by a microcanonical ensem- 
ble. When thermalized, the assymptotics of any expec- 
tation value of some operator A acting on the subsys- 
tem A can be obtained from statistical mechanical argu- 
ments. This is a topic of current intense research [l8|. 
A rule of thumb is that a chaotic system thermalizes in 
the above sense. The present model, in its closed coun- 
terpart, Eq. ([2]), is not chaotic and contains only a few 
degrees of freedom, and one thereby do not expect ther- 
malization. In fact, we know that for relatively localized 
initial states ((An) ~ (n)), the JC model possesses a se- 
ries of collapse-revivals and we will not have any 
thermalization. Things may become more interesting if 
we go beyond the RWA and include the counter rotating 
terms to our Hamiltonian. The model is then believed to 
be non-integrable [l3, E^l- Even so, within the range of 
our numerics, generic thermalization is not found in our 
model. Instead, we will consider an even weaker type 
of thermalization that was recently discussed by Pono- 
marev et al. in Ref. jsij . Their model does not put 
any constrains on integrability nor on a large number of 
degrees-of-freedom. In this respect it is more related to 
the system at hand. 
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FIG. 8: Evolution of the scaled photon inversion Z{i) (a) and 
the Bloch vector amplitude |-R(f)| (b) for an initial thermal 
state in mode a and vacuum in b. The inversion shows a 
Gaussian collapse, and at later times revivals build up indi- 
cating absence of thermalization. The Bloch vector amplitude 
collapses and revive at half the inversion revival and also at 
the revival of the inversion. In between these revivals a series 
of fractional revivals are seen. The dimensionless parameters 
are cjc = 54.3(7, Q. = ISOg, and Ua ~ 5. 

The Hamiltonian considered in [3l| has the form 

Hp^HA+HB+H^nt, (28) 
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where Ha and Hb act in subsystems A and B re- 
spectively, and they are assumed identical Ha — Hb, 
and furthermore the interaction is taken on the form 
Hint — Ya Yb where Ya — Yb and they too act 
in the respective subsystems. Thus far, the model is 
rather generic. However, one more constrain is assumed, 
letting En be the n'th eigenvalue of Hp they assumed 
En — Em = Eg — Ew if and only \i s = n and w = m. 
Within this model, Ponomarev et al. showed that for 
an initial state p(0) = pa(0) (8i Pb(0), with /3a,b(0) ther- 
mal states characterized by the temperature T4 and Tg 
respectively, the asymptotic states pa(oo) = P_b(oo). 
Hence, regardless of the two temperatures Ta and Tb, 
the system thermalizes to states with a common temper- 
ature. It was further demonstrated that this property 
is not restricted to initial thermal states, but holds for 
more general diagonal initial states pa{0) and /9b(0). 

As a closed system ([25)) . we will in this section neglect 
any irreversible coupling to reservoirs. Thinking of sub- 
system A and B as our two boson modes, it is clear that 
our model does not fulfill the above constrains. First 
of all, the qubit adds extra degrees-of-freedom, and sec- 
ondly, not even in the linear regime (i.e. the importance 
of the qubit degrees-of-freedom has been omitted) does 
our Hamiltonian possesses the desired form (^5)) . Fur- 
thermore, within the RWA, the linearized model has a 
spectrum that does not obey En — Em ^ Eg — E^ for 
different (n,m) and {s,w). This condition can, however, 
be met by including the counter rotating terms, in other 
words go beyond the RWA. To have a time-independent 
model we do not consider external pumping, and our 
Hamiltonian to be studied in this section is then 



H = LJc (ji^ a + b^b^ + i^^z 

+g (flt + a) + (h^ + bj a 



(29) 



We know that the qubit induces an effective coupling 
of the two modes, and that the counter rotating terms 
brings about an irregularity of the spectrum. Therefore, 
it feels natural to expect our non-RWA system to ther- 
malize in the sense of Ref. [SlJ . 

As our initial state we let the qubit be in the same 
pure state as earlier, mode b in vacuum, and mode a in 
a thermal state 



PA{t - 0) 



^ (n, + l)"+i 

n— 



\n) (7i|, 



(30) 



where determines the number of thermal photons, i.e. 
(?^a) — At zero temperature, T = 0, we naturally 
have fia — 0. Note that beyond the RWA, the number of 
excitations is not conserved. If we assume that the energy 
exchange between the fields and the qubit is small within 
the large detuning regime, thermalization then implies 



with n = fiajl- 

Our numerical results, presented in Fig. |S1 demon- 
strates a Gaussian collapse of the photon inversion into 
a state with a thermal distribution. However, after some 
time the Josephson oscillations sets off again and the sys- 
tem revive. The two collapsed field states showed identi- 
cal thermal like photon distributions, but the off diagonal 
coherence terms of the states are not exactly zero which 
causes the revival. In the JC model, the revivals become 
more fuzzy with time and after very long times the re- 
vivals begin to overlap causing super-revivals. The well 
resolved collapse-revival structure of Fig. |S] (a) seems to 
survive for extremely long times, and within our numeri- 
cal tests we have surprisingly not seen any long time de- 
cay of the revivals. The amplitude of the revivals varies, 
but every now and then a revival appears where Z(i) os- 
cillates approximately between -1 and 1. Furthermore, 
we have checked our findings for different initial popu- 
lations fia and A. However, if A is taken too small the 
non-linearity is very strong and the whole collapse-revival 
structure smears out. The Bloch vector amplitude dis- 
played in Fig. [5] (b) reproduces a similar collapse-revival 
structure as was found in Fig. [5] (b), i.e. consisting of 
a series of full and fractional revivals. Thus, once more 
does the qubit possesses a more complex dynamics than 
the Josephson current. The increased qubit purity at 
half the revival time is reminiscent of that encountered 
in the regular JC model Not shown in Fig. |8] is 

the field purity. Despite the fact that the initial state 
is mixed and one would expect thermalization, we have 
found that the field purity displays a similar structure as 
the Bloch vector amplitude. Thus, the structure of the 
field purity contradicts as well thermalization. We should 
point out that we have also verified similar results as the 
ones above when the qubit is an initial maximally mixed 
state, p,t{t - 0) = (|+)(+| + |-)(-|)/2. 

Figure [S] hints that our system does not thermalize. 
To explore this further, we consider instead an initial 
diagonal state 



PA 



\n){n\ 



(32) 



■11=0 



Thus, despite being diagonal the photon distribution is 
Poissonian as for a coherent state. In Fig. [5] we display 
the corresponding field linear entropy Pf(t) (a). This 
is compared with the field linear entropy with an initial 
pure state pA{t = 0) = |Q;)(a|. The general structure is 
the same between the two cases, but it is clear that even 
at the well resolved revivals the purity of the initially 
mixed state is not increased considerably. Nonetheless, 
it clearly demonstrate the absence of thermalization. 

Finally we explore a different definition of thermaliza- 
tion which was discussed in Ref. [13] ■ With the model 
Hamiltonian d^H]) , the time-averaged reduced density op- 
erators 



Pa(oo) ^Pb(oo) (31) 



— lim — 



dtpa (t), 



A, B, 



(33) 
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FIG. 9: Evolution of the field linear entropy Pf{t) for an 
initial mixed field state of mode a H32|) (a), and a coherent 
state with the same photon distribution as the mixed state 
(b). The 6- mode is initially in vacuum and the qubit as in 
previous examples. The general structure is similar in both 
cases, but the initial mixed state stays largely mixed even at 
the revivals in contrast to the initial pure state which almost 
perfectly regains its purity at the half and full revivals. The 
dimensionless parameters are as in Fig. [S] except the initial 
mean number of photons fia ~ 16. 



where analyzed. If pa becomes diagonal, the system 
clearly has thermalized. We do not show our results here, 
but we have verified for different initial states and param- 
eters that the time-averaged reduced density operators do 
not diagonalize. Thus, once again we have found absence 
of thermalization. 



VII. CONCLUSIONS 

We have presented a most simple model that simulates 
the photonic counterpart of a bosonic Josephson junc- 
tion. Opposed to earlier studies of the photonic Joseph- 
son effect, our model describes the internal effect. In 
the dispersive regime a qubit mediates an effective cou- 
pling between two photon modes of either an optical or 
a transmission line resonator. After introducing the gen- 
eral model we demonstrated the Josephson effect both 
within a linear analytic model and using full numerical di- 
agonalization. Our numerical results shows the presence 
of collapse-revivals originating from the intrinsic qubit- 
field non-linearity. 

We continued by addressing several aspects that arises 
from decoherence in one or another way. To start with, 
we considered inherent decoherence of the qubit stem- 
ming from interaction with the quantized boson modes. 
In such a Markovian analysis we found a rapid Gaus- 
sian decay of the qubit purity, at quarter of a Josephson 
oscillation period the qubit has lost almost entirely its 
coherence. At half an oscillation, the coherence is essen- 
tially completely regained. Our setup serves as well as a 
realization of the "quantum gatekeeper model discussed 
in (so]- As in Ref. [SOj, regarding a qubit coupled to a 



boson Josephson junction formed in a Bose-Einstein con- 
densate |30|] , we also examine the prospects of using the 
qubit as a quantum meter for the Josephson oscillations. 
We show that the Josephson current can be measured by 
integrating the length of the Bloch vector. 

What was not put forward in [30,] is the backaction 
a qubit measurement has on the Josephson current it- 
self. As a first example, we mimicked the irreversible 
qubit measurement within a master equation approach. 
The induced dephasing of the qubit rapidly propagates 
to the boson modes which leads to an exponential de- 
cay of the oscillations. The more information subtracted 
from the qubit, the stronger is the demolishing effect on 
the Josephson oscillations. This fact demonstrates the 
macroscopic nature of the Josephson oscillations as a col- 
lective quantum many-body phenomenon. As a second 
example, we considered non-selective or projective mea- 
surements, and such interrupted unitary evolution was 
shown to render a quantum Zeno slowing down of the 
Josephson oscillations. 

The following section looked into effects deriving from 
cavity dissipation. Again, this was modeled utilizing a 
master equation. The important result coming out from 
this analysis was the fact that even when realistic exper- 
imental parameters was considered, the Josephson oscil- 
lations survived for around 10 full oscillations. This is 
indeed comparable, or even longer, than what has been 
observed in either Bose-Einstein condensates trapped in 
double- well potentials [1^ , or i n sp inor condensates with 
Raman coupled Zeeman levels [61] . A somewhat surpris- 
ing finding was that photon dissipation did not render an 
exponential decay of the scaled Josephson current, but 
instead we found a Gaussian decay. This is in contrast 
to the exponential decay of the amplitudes of the two 
modes, and also to the exponential decay of the scaled 
current when exposed to dephasing. 

We finished our investigations by considering thermal- 
ization in a weak sense. That is, we did not study equili- 
bration into microcanonical ensembles where the expec- 
tation values (A) can be predicted by the correspond- 
ing ensembles [1^ , but thermalization of two interacting 
thermal quantum systems according to [sij . Our model 
fulfills most of the constrains for such thermalization, 
but it was demonstrated that despite this we do not find 
thermalization. Lack of the sort of thermalization as dis- 
cussed in [s^ l was also numerically verified. 

Throughout we assumed the two boson modes to be de- 
generate. Our results would still be valid even if we drop 
this constrain. Then, however, should it be noted that 
the Josephson effect is detuned and hence not all pho- 
tons will be swapped between the two modes. In circuit 
QED, which utilizes transmission line resonators, there is 
in general only a single mode for a given frequency. How- 
ever, in optical Fabry-Perot cavities there are typically 
two orthogonally polarized degenerate modes. Notewor- 
thy is that bimodal cavity QED experiments have already 
been demonstrated |63|]. 
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